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We consider nonintersecting random walks satisfying the condi- 
tion that the increments have a finite moment generating function. 
We prove that in a certain limiting regime where the number of walks 
and the number of time steps grow to infinity, several limiting dis- 
Qh ' tributions of the walks at the mid-time behave as the eigenvalues of 

pH , random Hermitian matrices as the dimension of the matrices grows 

* ■ to infinity. 

■ 

1. Introduction. It is known that various limiting local statistics arising 
in random matrix theory are independent of the precise structure of the 
J> ■ randomness of the ensemble [11, 17, 18, 21, 32, 44, 48]. For example, consider 

CN . the set of Hermitian matrices equipped with a probability measure invariant 

under unitary conjugation. For a very general class of measures, as the size 
{/-^ ■ of the matrix becomes large, the largest eigenvalue converges in distribution 

to the Tracy-Widom distribution, while the gap probability in the "bulk 
scaling limit" converges to a (different) universal distribution. 

It has been discovered that the limiting distributions arising in ran- 
dom matrix theory also describe limit laws of a number of specific mod- 
els in combinatorics, probability theory and statistical physics; apparently, 
these models are not expressible in terms of random matrix ensembles. 
Examples include the longest increasing subsequence of random permu- 
tations [5, 13, 31, 43], random Aztec and Hexagon tiling models [9, 33], 
■ last passage percolation models with geometric and exponential random 
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variables [30], polynuclear growth models [34, 45] and vicious walker mod- 
els [3, 27]. For these models, the distribution function of interest was com- 
puted explicitly in terms of certain determinantal formulae and the asymp- 
totic analysis of these determinants yielded the desired limit law. Neverthe- 
less, it is believed that such limit laws should hold for a class of models 
much wider than the explicitly computable ( "integrable" ) models. One such 
universality result for models "outside random matrices" was obtained in 
[10, 12, 50] for thin last passage percolation models with general random 
variables. 

This paper studies nonintersecting random walks and proves random ma- 
trix central limit theorems in a certain limiting regime. The motivation for 
this study comes from two sources. The first is the fact that the eigenvalue 
density function of the Gaussian unitary ensemble can be described in terms 
of a nonintersecting Brownian bridge process [22, 33]. Namely, consider n 

standard Brownian bridge processes (B^\ . . . , .B t ) conditioned not to in- 
tersect during the time interval (0, 2) (i.e., b\ 1 ^ > • • • > for < t < 2), 
all starting from and ending at the origin. A simple computation shows that 
the distribution of {B^\ . . . ,B^} at time 1 is the same as the distribution 
of the eigenvalues of the n x n Gaussian unitary ensemble; see Section 1.1.1 
below for the computation. Hence, it is natural to ask if the same limit laws 
hold for general nonintersecting random walks. The second motivation is 
that a number of the aforementioned specific probability models for which 
the random matrix central limit theorem was obtained are indeed inter- 
preted in terms of nonintersecting random walks. We mention a few of them 
in the following subsection. 

1.1. Motivating examples. We begin by introducing two distribution func- 
tions. Define the kernels 

a — b ir(a — b) 

Set 

(2) fJ rw (0=det(l-A| (e>oo) ), F Sin e(r/) = det(l-S| [ _ r?j7?] ). 

The Tracy-Widom distribution, -Ftw, is the limiting distribution of the 
largest eigenvalue and Fsine is the limiting distribution for the gap proba- 
bility of the eigenvalues "in bulk" in Hermitian random matrix theory. 

1.1.1. Nonintersecting Brownian bridge process. Let Bt = (B^ , . . . , B^) 
be an n-dimensional standard Brownian motion. We compute the density 
function of B 1 conditioned that B^ > B ( t 2) > • • • > B { t n) for < t < 2 and 
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B =B 2 = (0,.. . ,0). Let pt(x,y) = ^ t e ~ {x ~ y)2/{2t) ■ The argument of Kar- 
lin and McGregor [36] implies that the density function of n one-dimensional 
nonintersecting Brownian motions at time t which start from (sci, . . . ,x n ), 
where x\ > ■ ■ ■ > x n , is given by 

(3) f t (b 1 ,...,b n ) = det(p t (xi,bj))i tj=1 , b 1 >--->b n . 

Hence, for b\ > ■ ■ ■ > b n , the density function of B\ equals 



(4) 



f(b 1 ,...,b n )= lim 

x,y-*o aet{p2{x i ,y j ))fj =1 

2»i(r»-l)/2 n 
= Wmn-i, II l^-M 2 IK 

11 3 =1 J- l<j<,< n i-l 



l<j<j'<n j — 1 

which is the density function of the eigenvalues ofnxn Hermitian matrices 
from the Gaussian unitary ensemble. Therefore, combined with the well- 
known results of random matrix theory, 

(5) lim F((b\ 1] - V2^)V2n 1/6 <x) = F TW (x). 

n — >oo 

1.1.2. Longest increasing subsequence and Plancherel measure on parti- 
tions. The longest increasing subsequence problem can be formulated in 
the following manner. Denote by S n the symmetric group on n symbols en- 
dowed with uniform measure. Given ir £ S n , a subsequence n(h), ■ ■ ■ , 7r(v) 
is called an increasing subsequence if %\ < ■ ■ ■ < i r and 7r(ii) < • • • < 7r(v). 
Denote by l n {^) the length of the longest increasing subsequence (this sub- 
sequence need not be unique). For applications of l Tl and activities around 
the asymptotic behavior of £ n , see, for example, [2, 5, 16]. In particular, 
in [5], the following limit theorem is: 

(6) lim pf W ~ 2 ^ <x)= Frw(x). 

A closely related object is the uniform measure on the set of pairs of 
standard Young tableaux having the same shape (equivalently, the so-called 
Plancherel measure on the set of partitions). Given a partition of n, A = 

(Ai, . . . , A r ), where Ai > • • • > A r > and Ai H h A r = n, a standard Young 

tableaux of shape A consists of r rows of boxes with distinct entries from 
{1, . . . ,n} such that the rows are left-justified, the ith row has Aj boxes and 
the entries are constrained to increase along rows and columns from left 
to right and top to bottom, respectively. These objects will be called row- 
increasing Young tableaux if the rows increase but the columns do not neces- 
sarily increase. The Robinson-Schensted bijection implies that the number 
of boxes in the top row of the pair of standard Young tableaux correspond- 
ing to 7T G S n is equal to £ n (ir) [49]. Therefore, the distribution of £ n is 
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the same as the distribution of the number of boxes in the top row of the 
pair of standard Young tableaux having the same, shape chosen uniformly. 
This correspondence provides a representation of l n which is computable in 
terms of explicit formulae if the number of standard Young tableaux of a 
given shape is computable. 

One way (among many) to compute the number of standard Young tableaux 
of shape A is by means of a nonintersecting path argument [35] . Let A^ 1 , . . . , N£ 
be independent rate-1 Poisson processes with initial conditions iVg = 1 — i 
for % = 1, 2, . . . , r. Define A\ to be the event that N[ = A.; + (1 — i) for all 
i = 1,2, ...,r. For almost every element of A\ (the elements of A\ where 
no two jumps of these processes occur at the same time), there is a natural 
map to a row-increasing Young tableaux. The map is defined as follows. If 
N l jumps first, then place a 1 in the leftmost box in row i; if N J jumps 
second, then place a 2 in the first box of row j if j ^ i and a 2 in the second 
box of row i if j = i. Continue in this fashion to produce a row-increasing 
Young tableaux of shape A. It is not hard to show that this map induces the 
uniform probability measure [when properly normalized by P(^4a)] on the 
row-increasing Young tableaux. The subset B\ C A\ which is mapped to the 
standard Young tableaux of shape A corresponds to the set of realizations 
whose paths do not intersect each other for all t £ [0, 1] . Since the mapping 
described induces uniform measure on the row-increasing Young tableaux 
of shape A and the standard Young tableaux correspond to nonintersecting 
path realizations, B\, the number of standard Young tableaux of shape A, 
can be computed by evaluating 

(7) | row-increasing Young tableaux of shape ^I p7^ y 

The denominator of (7) is e~ r Y\l =1 pj, by definition of Poisson processes 

and the independence of the N l , while | row-increasing Young tableaux of 
shape A| = A , by elementary combinatorics. On the other hand, via the 
Karlin-McGregor formula [36], 



(8) p ^)= de \ (A ,_ i+i) , 

Hence, the number of standard Young tableaux of shape A is 
n ' det( a )i,j=i ■ I n tandem with the RSK correspondence, this formula 

leads to an algebraic formula for the number of ir G S n for which £ n (Tr) < m. 
Moreover, a slight extension of this argument shows that result (6) can be 
stated in terms of the top curve of the nonintersecting Poisson processes if 
these processes were forced to return to their initial locations at time 2 by 
imposing that their dynamics between times 1 and 2 have negative rather 
than positive jumps. The asymptotic behavior of other curves can also be 
studied [6, 7, 13, 31, 43]. 
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1.1.3. Symmetric simple random walks and random rhombus tilings of a 
hexagon. Consider n symmetric simple (Bernoulli) random walks S(m) = 
(S^ (m), ■ ■ ■ , S^ n '(m)), conditioned not to intersect and such that S(0) = 
(2(n — 1), 2(n — 2), . . . , 0) = S(2k). Any realization of such walks is in one- 
to-one correspondence with a rhombus tiling of a hexagon with side lengths 
k,k,n,k,k,n. Again, using the argument of Karlin and McGregor, the distri- 
bution of S(k) can be expressed in terms of a determinant. This determinant 
was significantly simplified and was shown to be related to the so-called Hahn 
orthogonal polynomials by Johansson [33]. A further asymptotic analysis of 
the Hahn polynomials [8, 9] shows that as n, k — > oo such that k = 0(n), 
the top walk S^\k) converges to Ftw and the gap distribution "in bulk" 
converges to a discrete version of i*Sine • A similar asymptotic result was also 
obtained for domino tilings of an Aztec diamond [33]. 

Certain polynuclear growth models, last passage percolation problems 
and a bus system problem [4, 34, 42, 46] have also been analyzed in depth 
using nonintersecting path techniques. In each of the cases described above, 
the random walks are very specific and the analysis relies heavily on their 
particular properties. 

1.2. Statement of theorems. Let k be a positive integer. Let 

(9) Xi = ^-=A ie{0,. ..,*:}. 

Note that Xi 6 [— 1,1] for all i. Let {i^ J }^lf^^ =1 be a family of indepen- 
dent identically distributed random variables where is a positive integer. 
Assume that Wfj 1 = and Var(Y^) = 1. Further, assume that there exists 

A > such that E(e XY i J ) < oo for all |A| < A . 

Define the random walk process S(t) = (So(t), . . . , Sk(t)) by 



[Tf lt ^l 21 i ftN k tN k 

do) *,-(*) =* i + ^ p + 

for t £ [0,2], 

which starts at Sj(0) = Xj. For N k equally spaced times, Sj is given by 

(U) Sj ilk 1 ) = Xj + {^k {Yl + " ' + Y ' )j 1 = 1 > 2 >-> N »- 

For t between -j^-l and + 1), Sj(t) is simply defined by linear interpo- 
lation. 

Let (C([0,2];R fc+1 ),C) be the family of measurable spaces constructed 
from the continuous functions on [0, 2] taking values in equipped with 
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their Borel sigma algebras (generated by the sup norm). Let A k ,B k £C be 
the events defined by 

(12) A k = {y (t)<---<y k (t) for t€ [0,2]}, 

(13) B k = {yi{2) £ [xi - h k ,Xi + h k ] for i G {0, . . .,k}}, 

where h k > 0. The results of this paper focus on the process S(t) conditioned 
on the event A k C\ B k , where <S f ■ In other words, the particles never 
intersect and all particles essentially return to their original locations at the 
final time 2. 

The main results of this paper state that under certain technical condi- 
tions on h k and N k , as k — > oo, the locations of the particles at the half time 
(t = 1) behave statistically, after suitable scaling, like the eigenvalues of a 
large random Hermitian matrix from the Gaussian unitary ensemble. The 
conditions for h k and Nk are that {h k } k> o is a sequence of positive numbers 
and that {N k } k> o is a sequence of positive integers satisfying 

(14) h k <(2k)~ 2k2 and N k >h~ A(k+2) . 
Let C k , D k £ C be defined by 

(15) C k = {y k (l)<V2k + 7 ^], 

(16) £>* = {l/i(l)£ [-^,^] foralHG{0,...,fc}|, 

where ^ and 77 > are fixed real numbers. The event C k is a constraint on 
the location of the rightmost particle and D k is the event that no particle 
is in a small neighborhood of the origin at time 1. 

Theorem 1 (Edge). Let F k be the probability measure induced on (C([0, 2]; 
M fc+1 ),C) by the random walks {S(t):te [0,2]}. Let {h k } k>0 and {N k } k>0 
satisfy (14)- Then 

(17) lim F k (C k \A k n B k ) = F TW (0- 

k— >oo 

A similar theorem holds for the bulk. 

Theorem 2 (Bulk). Let ¥ k be the probability measure induced on (C([0, 2]; 
M fc+1 ),C) by the random walks {S(t) :t G [0,2]} and let {h k } k> o and {N k } k> o 
satisfy (14)- Then 

(18) lim P k {D k \A k n B k ) = Fs^). 

k — >oo 
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The proofs have a two-step strategy. The first step is to show that under 
the conditions of the theorems, the process S(t) is well approximated by 
nonintersecting Brownian bridge processes starting and ending at the same 
positions. This proof relies on the Komlos-Major-Tusnady (KMT) theorem. 
The second step is to compute the limiting distributions of the noninter- 
secting Brownian bridge processes and prove that these distributions are 
indeed i*rw or Fsi ne . This process is quite similar to the one discussed in 
Section 1.1.1, with the minor change that the Brownian bridge processes 
start and end at equally spaced locations, rather than at the same location. 
This change results in a Coulomb-gas density with the so-called Stieltjes- 
Wigert potential instead of the quadratic potential which appears in the 
GUE case. Such a nonintersecting Brownian bridge process was also con- 
sidered in [24, 26] and the connection to the Stieltjes-Wigert potential was 
made in [26] in order to compute the partition function and the limiting 
density of states. However, the edge and bulk scaling limits of the system 
had not been worked out. This paper obtains the asymptotics of the or- 
thogonal polynomials with respect to the Stieltjes-Wigert weight by using 
the Riemann-Hilbert method. As a consequence, the edge and bulk scaling 
limits are obtained. 

The above theorems are proved under the condition that N). is large 
compared to k + 1, the number of particles. This assumption ensures that 
the Brownian approximation of the random walks has a smaller effect than 
the nonintersecting condition. Although it is believed that the condition 
on Nk is technical, it is not clear under which conditions on the random 
variables one has N). = 0(k). For example, when {Y^} are Bernoulli, these 
results were proven even when Nk = 0{k) (see Section 1.1.3 above). This 
is because there is an integrability in this problem: the Karlin-McGregor 
argument applies directly because intersecting paths must be incident at 
some time. It is a challenge to find the optimal scaling such that a result 
of this nature holds for more general random variables. In other words, in 
what scaling regime does the exact Karlin-McGregor calculation essentially 
not matter? 

This paper is organized as follows. The approximation by the Brownian 
bridge process is proved in Section 2. The asymptotic analysis of the Brow- 
nian bridge process (appearing in Section 2) is carried out in Section 3. 
Some other considerations such as finite-dimensional distributions and the 
modifications necessary to study random variables without finite moment 
generating functions are discussed in Section 4. 

2. Approximation by a Brownian bridge process. Let {Xt}t>o be the 
M fc+1 -valued stochastic process X t = (Xo(t), . . . , Xk(t)), where Xj(t) = Xj + 
B{ for a family of k + 1 independent standard Brownian motions B\ . The 
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proof in this section relies on the Komlos-Major-Tusnady coupling of Brow- 
nian motions and random walks [38, 39] which can be stated in our set- 
ting as follows. With increments of the form {Y^} j' = q 1=1 described in the 
Introduction, there exists a coupling such that 



(19) P 



sup 

0<l<N h 



Si 



21 



Xi 



21 



> 



1 



[c\ogN k + x) <e 



for some fixed a, c > which depend only on the properties of the moment 
generating functions of the {Y^} k -^Q l=l . Alternatively, (19) can be written 

as 



(20) 



sup 

0<l<N k 



Si 



21 



Xi 



21 
N~k 



dogN k 
> — 7=^ + y 



< e 



This fact immediately implies that 



(21) 



sup sup 

0<i<kO<l<N k 

<(k + l)e- ay ^ 



Xil — 



2/ 



clogN k 

> — t=- + y 



Let {<5(i)}te[o,2] be the (k + l)-dimensional random walk process defined 
in the Introduction and let {A^} be the KMT-coupled (k + l)-dimensional 
Brownian process on the same probability spaces 

(ft( fc ) | ,F(*),p( fc )). We can 

assume that the probability space which holds S and X is large enough to 
hold a third process Z t = (Zo(t), . . . , Z k (t)), where the Zi(t) are standard 
Brownian bridge processes with initial and terminal conditions specified by 
Zi(0) = Zi{2) = Xi . Let F fc 5 , F x , Fj? : C([0, 2], M fc+1 ) ^1 be defined by 

^-A k nB k nC k (y) 



(22) 
(23) 



F k x (y) 



F k z (y) 



(24) 

Let Gf,G^,Gf:C([0,2],R fc + 1 )- 

(25) G s k (y) 

(26) G${y) 

(27) G z k {y) 



El AknBk (S) ' 

^-A k nB k nc k (y) 
Et AknBk (X) ' 

^A k nB k nc k {y) 

m AknBk (z) ■ 

1 be defined by 

^-A k nB k nD k (y) 
VtA k nB k (S) ' 

^-A k nB k nD k (y) 
Et Ak nB k (X) ' 

^A k nB k nD k (y) 
Et Ak nB k (Z) ' 
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Theorems 1 and 2 will be proven in two steps. The first step is to show 
that under the conditions given in the Introduction, the following holds. 

Proposition 1. As k — ► oo, the random variables Fj? , Fj 2 , Gf and G z 
satisfy 

(28) E(if(S)-if(Z))-0, 

(29) E(G s k (S)-G z k (Z))^0. 

As E(if (5)) = P fe (Cfcl^ n fl fc ) and E(G?f(5)) = P fc (£> fc |A fc n 5,.), it is 
enough to prove the following. 



Proposition 2. 
(30) limEif (Z)=F TW (£), 



f( 



(31) lim EGt(Z)=F Sine ( V )- 

fe^oo 

The proof of Proposition 2 is given in Section 3 below. The rest of this 
section focuses on the proof of Proposition 1. Proposition 1 is proved in two 
steps: first, E(Ff (5)) is approximated by E(Fj*(X)) and second, E(F fc x (X)) 
is approximated by E,(F k (Z)). The proof of (29) is handled in a similar way. 

Three preliminary lemmas are needed in order to prove Proposition 1. 
Recall from (9) that 

(32) Xi = 2% ^ k , i = 0,...,k. 

Lemma 1. Let a, c > be the constants in the KMT approximation (21 ). 
For any p > 3cl ^ , 

E\l AknBk (S) - t Ak nB k (X)\ < 2(fc + l) e -(V2)aV^P 
(33) 

+ Me-W" + 8(2* + l)p, 

P V T 



(34) 



(35) 



E|lB fc (S f ) - ±B k (X)\ < (k + l)e- {1/2)aVNkp 

16(fc + l) /JV*. 



+ 



E|l Cfc (5) - l Cfc POI < (fc + l) e -(V2)«^P 



^K e -P 2 N k M + 8{k + l)()j 



+ 16(fc + l) /iV, e _ p2jVfc/64 + 
P V vr 
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(36) 



E\l Dk (S)-l Dk (X)\<(k + l)e 



-(l/2)a^N k p 



+ 



P V 7T 



(37) 



Proof. Note that 

E\± AknBh (S)-l AhnBk (X)\ 

= E\l Ak (S)l Bk (S)-l AK (X)l Bk (X)\ 

= E\(l Ak (S) - l Ak (X))t Bk (S) + (l Bh (S) - l Bk (X))l Ak (X)\ 
<E\l Ak (5) - l Ak (X) \+E\t Bk (5) - l Bk (X) | . 



We first estimate E\t Ak (S) - lA k (X)\ = F{£), where £ = {S G A k , X <£ A k }U 
{S i A k , X G A k }. Recall that A k = {y (t) < ■ ■ ■ < y k (t) for t G [0,2]}. Let 



P> 



3clogjV fc 



, where c is the KMT coupling constant. The event £ can be ex- 



pressed as the disjoint union of the three events £1,82, £3- The first event £\ is 
the subset of £ consisting of "bad" paths satisfying sup 0< j <fc sup tg [ 0j 2] \Si(t) — 
Xi(t)\ > cl ^^ fc + p. The second event £2 is the subset of £ \ £\ consisting 

of paths satisfying min tg r 0] 2](5'j(t) — 5j_i(t)) < for some 1 < i < k while 
X (t) < ■■■ < X k (t) for all t G [0,2]. The third event £ 3 is the subset of 
£\£\ consisting of paths such that min tg j 0i 2] ~~ -^i-i(i)) < for some 

1 < i < k while S (t) < ■ ■ ■< S k (t) for all t g'[0, 2]. In order to estimate P(£i), 
note that the KMT theorem couples random walks to Brownian motion at 
discrete times. Hence, even when X and S are close at discrete times, "bad 
paths" may occur if X fluctuates too much in (j^l, jr{l + 1)) for some I. 
(Note that S is simply linearly interpolated for times not integral multiple 
of jj--) Thus, from (21) and standard estimate for Brownian motions, 

P(£i)<p(su P sup Si(^P)-X 



sup sup 

0<i<kO<l<N k 



+ 



sup sup 

0<i<kO<l<N k 



Si 



21 



21 \ 
N~k) 

X, 



dogN k p 

2 



21 



n 



max \Xi(t) 

s,t&(2l/N k ,2(l+l)/N k ) 



< clogiVfc p 

P 



(38) 



for some < i < k and for some < I < Ni 



< (k + l)e 



-(l/2)a^N k -p 



+ (k + l)N k F[ max {XUt) - XUsM > £ 

\t,s£[0,2/N k ] 2 
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(l/2)av^p , 16 ( fc + 1 ) t / -fP Nk /U 



P V I" 

Note that this estimate does not use the fact that £\ is a subset of £. For a 
path in the event £2, there exists z € {1, 2, . . . , k} such that min tg r 0] 2](5'j(t) — 
Si-i(i)) < 0, but < for all t G [0,2] and |S,-(i) - < 

Cl /Af Vfc + P all i G [0, 2] and j € {0, 1, ... , ^Vfe}- Therefore, for a path in £2, 

v ^ k 

< mm te[0i2 ](Xj(t) - AVi(t)) < 2c ^ ffc + 2p < 4p. Thus, from a standard 
Brownian motion argument, 



2) < P < min (XAt) - Xi-At)) < 4p for some 1 < i < k 
\ te[o,2] 
(39) 

< fcPj < min (Xi(t) - X (t)) < 4p] < 4k p. 
V te[o,2] / 

A similar argument yields that 



(40) 



3 j < P -4/0 < min - AVi(i)) < for some l<i<k 

te[o,2] 



< MP -4/0 < min (Xi(t) - X (t)) < < 4kp. 
te[o,2] 



Therefore, 

E|l Afc (£)-l Afc (X)| 
(41) =P(£ 1 )+P(£ 2 )+n 



< (fc + i) e -(V2)«^P + 16(fc + l) Nk e - p iN k /u + 

p V 7T 

We now estimate E|l Bfc (5) - = P(J zr ), where f = {5eB fc ,I^ 

.B/t} U {S ^ B k ,X G -Bfc}. As before, we express T as .Fi U T2 U JT 3 , a disjoint 
union. The first event, Fi , is the subset of T consisting of the same bad paths 
as in £\. The event T2 is the intersection of T \ T\ and {S G B k , X ^ B k } 
and the event J-3 is the intersection of J-\J-\ and {S ^ B k ,X G -B^}. The 
argument for £\ implies that the same bound (38) applies to P(.Fi). For a 
path in T2, there exists i G {0, 1, . . . , k} such that Xi(2) ^ [scj — + fofc]. 
But as 5,(2) G [2* - h k , Xi + fc fc ] and \Si(2) - A,(2)| < ^£|k +p <2p, we 

find that X;(2) G (xi + h k ,Xi + h k + 2p] or Xi(2) G [xi - h k - 2p,Xi - h k ). 
Therefore, 

(42) P(F 2 ) < (fc + 1)P(A (2) G [-2p, 2p]) < 4(* + l)p. 

A similar argument yields the same bound for ¥(J-s). Therefore, (34) is 
proved, as is (33), by using (37) and (41). An almost identical argument 
proves (35) and (36). □ 
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Denote by pt{a,b) = - 7 =e~( a_ '') 2 /( 2t ) the standard heat kernel in one 
dimension. The theorem of Karlin and McGregor [36] for nonintersecting 
Brownian motions implies that the joint probability density function 
ft(yo, ■ ■ ■ ,Uk) of {k + l)-dimensional Brownian motion X(t) at time t sat- 
isfying Xq(s) < Xi(s) < • • • < Xk(s) for s £ [0, t] is equal to 

(43) f t (y , ...,y k ) = det(p t (xi,yj)) k d=0 , 

where Xi = Xi(0). The following lemma establishes a lower bound for this 
density when yi = X{ for all i. 

Lemma 2. For t > 0, 

1 / 9 \ k(k+l) 

(44) det^,,,))^ > (^TW e - 2{fc+1)(fe+2)/(m ° (^) • 
In particular, for all sufficiently large k, 

(45) det(p 2 (x l ,x j ))l j=0 >k- k2 . 

Proof. As X{ = we have 

det^,^))^ = detf ie-VWW' 

(46) 



27rt / i j=o 

2 E^ 2 



det(e 



2ij/(tk 2 )\k 



(2iTt)( k + 1 )/ 2 """^ H ' j=0 ' 
It is an exercise to show that for k > 1, 



(47) det(e 2 ^/(*fc 2 )) fc 



i,3=0 



A: 



Jj5'(i-D/s 



A- 



J(,J3-1)*+1-J 

■i=i 



u=i 

where 5 = e 4 ^* fc2 \ Using (47) and the fact that 5 — 1 > vp- > 0, we obtain 



det(pt(xi,^))£ =0 



tfe 2 



1 k 

-2(fc+l)(fc+2)/(3tfc) TT^j _ ■QAs+l-j 



(48) 



(2vrt)( fc+1 )/ 2 



> 1 r -2(k+l)(k+2)/(3tk) ( r -i \fc(fc+l)/2 

" (2vrt)( fc + 1 )/ 2 1 j 

> 1 -2(fc+i)(fc+2)/(atfc) / 4 \ 
- (2vrt)( fc + 1 )/ 2 V** 2 / 



This completes the proof of Lemma 2. □ 
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The following lemma will be used to control the difference between a 
conditioned version of the process X and the process Z. 



Lemma 3. If h k < (2fc) 2fc2 ; then for sufficiently large k, 

rl f n-i f n- nr ■ 11 f 

(49) 



det(pi(7/i,xj)) " ' det(pi(yi,a; j + a,-)) ds - • ■ ds k 



det (pa (x< , Xj )) . . . j^ hk det (pa (a* , Xj + Sj )) ds - ■ ■ ds k 

uniformly in (yo, . . . , y k ) £ . 



1 



Proof. The conclusion of this lemma is a consequence of several ele- 
mentary determinant estimates. First, note that if A = ((%)f j=o is a (fc + 
1) x (k + 1) matrix with entries |aj,-| < 1, then for the matrix P 3 given by 
(I lj )mn = 5i m 5 jn , we have 



(50) 



|det(A) -det(A + eI ij )\ <ekl 



(51) 



Using a Lipschitz estimate for the Gaussian density, equation (50) implies 
that for any t > ^= , any /i > and any (ao, • ■ • , a k ), (&o 5 • ■ • ? fyfc) G M fc+1 , 

1 / >/l 

det(pi(Qj,6j)) - T^p+i y ft det (K( a i + s ii^)) rfs o---^s fe 

<2/i(/fc + l) 2 A;!. 

A simple algebraic manipulation now yields that the left-hand side of (49) 
equals 

(52) det (Pi(^' x j)) ^ 2 . Ql 

where 

(53) 



+ 



det(p2(xi,Xj)) det(p 2 (xi,Xj)) + Q 2 det(p 2 (xj, Zj)) + Q 2 



Si 



Q-2 



(2h h ) 



fc+i 



[-hfc,A*] fc+1 



(54) 



[det(pi(y;,x,,)) 

- det(pi(j/j,Xj + Sj))] ds ---ds k , 
[det(p 2 (xi,Xj + sy)) 



(2/i fc ) fc + 1 7[_ hfci ^ ]fc +i 

- det(p 2 (xj,Xj))] ds • • • ds k . 
Using the estimates (45) and (51), we obtain 
(55) det(p 2 (x i ,x J )) + Q 2 > k~ k ' 2 - 2h k (k + l) 2 k\ > ±£T fc2 

for sufficiently large k. Hence, again using (51), 



(56) 



Qi 



det(p 2 (xi,Xj)) + Q 2 



<2h k {k + lfk\k k . 
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On the other hand, as det(pi(xi,yj)) is the density function for (yo, . . . , y k ) E 
M> +1 , where M> +1 = {(yo, . . . , y k ) E R k+1 : yo < ■ ■ ■ < yk} corresponding to the 
probability of k + 1 Brownian motions starting from (xo, • ■ • ,x k ) and ending 
at (yo, . . . , yk) at time 1 without having intersected, it is clearly less than the 
same type of probability density function when a nonintersection condition 
is not imposed. Therefore, 



(57) 



and hence 



1 



&et(pi(xi,yj)) < Y[ —?==£ 

■ Z7T 



det(pi(yj,x,)) 



i=0 



Q-2 



-(l/2)(x,-y,) 2 < I 



<2h k (k + lfk\k 



2k 2 



(58) 

det(p 2 (xi,Xj)) det(p 2 (xi,Xj)) + Q 2 
Since /i^ is assumed to be less than or equal to (2k)~ 2k2 , (49) follows. □ 

PROOF of Proposition 1. Two estimates will be needed. Note that 

(59) \E(F k s (S)-F k z (Z))\<E\F k s (S)-F k x (X)\ + \E(F* (X) - if (Z))|. 

The first term on the right-hand side of (59) is estimated as follows: 

E\F^(S)-F k x (X)\ 

^-A k nB k nc k (S) t AknBknCk (X) 



E 



(60) 



ElA k nB k (S) El AknBk (X) 

= E\ (l Ak nB k nc k (S)(El AknBk (X) - El AknBk (S)) 
+ (^A k nB k nc k (S) - l Ak nB k nc k (X))El AknBk (S)) 

x (El AknBk (S)El AknBk (X))- l \ 



< 



< 



\El AknBk (X)-El AknBk (S)\ 
Et AknBk (X) 



+ E 



(&A k nB k nC k (S) - t Ak nB k nc k (X)) 



El AknBk (X) 



2E\l AknBk (S) - t AknBk (X) \+E\t Ck (S) - t Ck (X) | 
^U k nB k (X) 



By setting p = N, 1//4 in Lemma 1, for sufficiently large k, it is easy to 



check that 



(61) 



E\t AknBk (S)-t AknBk (X)\ 



< 



20k 



N. 



1/4' 



E\l Ck (S)-l Ck (X)\< 



k 
20k 



N. 



1/4' 
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On the other hand, by using (45) and the argument leading to (55), for 
sufficiently large k, 

Et Ak n Bk (X) 

= (2h k ) k+1 det(p 2 (x i ,x j )) 

+ (El AknBk (X) - {2h k ) k+l det{p 2 (x i) x J ))) 

= {2h k ) k+l det{p 2 (x i ,x j )) 



(62) 



+ 



-h k ,h k } k+1 



(det(p 2 (xi,Xj + Sj)) -det(p 2 (xi,Xj)))ds ■■■ds k 



> 



(2fr fc ) fc+1 
2k k2 



Hence, from (61), for sufficiently large k, 



(63) 



E\F£(g)-F?(X)\< 



120/c fc2+1 
(2h k ) k + 1 N 1 k /4 







as k — > co. For the second term of (59), note that the Karlin-McGregor for- 
mula for nonintersecting Brownian motions implies that [cf. (43) above] the 
density function of the nonintersecting Brownian bridge process Z evaluated 
at time 1 is equal to 

, CA s ,, \ det(p 1 (a^,y J 0)*jM) drt (Pi(w> a; j))*j==o 
(64) f(y , ...,y k ) = . 

Similarly, the density of the nonintersecting Brownian motion X evaluated 
at time t is equal to 



(65) 



f(yo,---,Vk) 



h-h k ,h k \ k+1 det ( pi ( Xi > yi))iJ=o det(pid/i, Xj + Sj)) k j=0 ds ■ ■ ■ ds k 



I[-h k ,h k ] k+1 det(P2(xi,Xj + Sj))ds - ■ -ds k 



Therefore, 

\E(F k x (X)-F k z (Z))\ 



E 



< 



(66) 



fc+i 



^-A k nB k nc k (X) _ t AknBknCk (Z) 
Et AknBk (X) ' Et AknBk (Z) 

det(pi(xi,yj)) 



l-h k ,h k ] k + 1 



x det(pi(yi,Xj + sj)) ds - ■ ■ ds k ^j 
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det(p2(xi,Xj + Sj)) dso--- ds k 

l[-h k ,h k \ k+1 



-1 



det {pi(xj,yj)) det (p 1 (y { , xj ) ) 



det(p 2 (xi,Xj)) 
By using Lemma 3 and (57), this implies that 



dyo---dy k . 



(67) 



\E(F*(X) - if (Z))\ < \ [ IdetCpiCx*,^))! 



dyo • ■ • dy k 



< 



k Jm k+1 



n 



-(xi-Vi) 2 /2 



> Li=0 



dyo • • • dy fc < 



The proof of (29) is exactly the same. This completes the proof of Proposi- 
tion 1. □ 



3. Asymptotics of a Brownian bridge process. We prove Proposition 2 
in this section. Together with the results of Section 2, this completes the 
proofs of Theorem 1 and Theorem 2. 

From the density formula of Karlin and McGregor for a nonintersecting 
Brownian bridge processes [36] [cf. (43)], 

E(if(Z)) 



(68) 



det(p2(si,Xj))f, =0 



i [det( Pl (x i ,y i ))t=o] 2 EEC 1 " Hiiytfdyj, 

j=0 



where Xi = and 

(69) 

Also, 



H 1 (y) = l 



,oo) 



(»)■ 



(70) 



E(Gf(Z)) = — — ^g— 



x / [det (p! ( s< , % ) )^ =0 ] 2 J] ( 1 - H 2 ( Vj ) ) d Vj , 

JR > j=0 



where 

( 71 ) w 2 (i/) = i^^/^fc+T^/^+TjCy)- 

We need the limit of (68) and (70) as fe — ► oo. 
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In the discussion below, TL{y) denotes either 7ii or Ti.2- Indeed, the alge- 
bra below works for arbitrary bounded functions 7i(y). Using the formula 
for p t and the definition of Xi . £111 elementary algebraic manipulation using 
Vandermode determinants yields that (68) and (70) are equal to 

k 

0<i<j<k j=0 



k 

(72) C' k .J h+i I] ^ Vllk -e^^ffil-Wye-^ d Vj 



where C k is the normalization constant so that (72) becomes 1 when Ti{y) 
0: 

p -(fc+l)(fc+2)/(3fc) 

(73) C' k 



det(p 2 (x i ,x j ))l j=0 (k + 1)1(2^+^ 

Note that the integration domain is changed to ]R fc+1 by using the symmetry 
of the integrand. A similar calculation for the case when the assumption 
that the particles start at equally-spaced locations and arrive at (yo, . . . ,y k ) 
at time 1 (without any assumption regarding what happens after time 1) 
can be found in equation (4.7) of [24]. Introducing the change of variables 
yj = -|log?Xj — 1, (72) becomes 

k 

^Cu-f J] (n,-^) 2 n(l-H(%))^-( fc2 / 4 )^^ 2 ^, 
JR + 0<i<j<k j=o u i 

where 

(75) W(«) = wQlog«-l 
and the normalization constant is given by 

£.fc+l„-(fc+l)(fc+2)/(3fc) 

(76) C k 



det(p 2 (x i ,x j ))l j=0 (k + l)l(27r) k + 1 ' 

This is the standard (5 = 2 ensemble in random matrix theory on the real 
half-line IR + with the weight 

(77) w ( u \ = I e -(fc 2 /4)(log«) 2 = e -(fc 2 /4)(log«) 2 -log«_ 

u 

Note that w(u) = o(u~ m ) for any m > as u — > +oo, and w(u) = o(u m ) for 
any m > as u { 0. 

With the change of variables u = e _2//fc x, (77) becomes 

(78) w{u)du = c-e-^'^°^ 2 dx, c = e 1 / fc2 . 

This is, up to a constant, the Stieltjes-Wigert weight, which is defined as 

(79) Tr-^e-^iogz) 2 
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(see, e.g., Section 2.7 of [51] or Section 3.27 of [37]). The moments for the 
Stieltjes-Wigert weight is an example of an indeterminate moment prob- 
lem; hence, there are several weights that have the same moments as the 
weight (79). Another interesting feature of the Stieltjes-Wigert weight (79) 
is that the corresponding orthogonal polynomials (called Stieltjes-Wigert 
polynomials) are examples of so-called (/-polynomials with q = e~ 1 ^ 2k ) (see, 
e.g., Section 3.27 of [37]). The connection between the nonintersecting Brow- 
nian bridge process Z and the Stieltjes-Wigert weight was first observed in 
[26]; the Stieltjes-Wigert weight also appears in [28]. 

Various (3 = 2 matrix ensembles of the form (74) (on both the real line 
and subsets of the real line) have been analyzed asymptotically and it has 
been proven that the local statistics of the "eigenvalues" (or the particles 
uq, . . . , Ufc) are generically independent of the potential w. For example, such 
"universality" is proved when w(x) = e~^ k+1 ^ v ^ for an analytic weight V 
on M. or M + satisfying certain growth conditions as x — ► =too (and as x — ► 
for weights on R + ) (e.g., [11, 21, 40, 44]) and when w{x) = e ~Q( x )^ where 
Q{x) is a polynomial (e.g., [20]). However, the asymptotic analysis of the 
ensemble with the weight given in (77) above does not seem to appear in 
the literature. It is well known [see (80) below] that the asymptotics of 
[3 = 2 ensembles amount to the asymptotic analysis of the corresponding 
orthogonal polynomials. For our case, we need the asymptotics of the or- 
thogonal polynomials of degree k and k + 1 with respect to the weight (77) 
as k — > oo; note that the weight also varies as k increases. The asymptotics 
of Stieltjes-Wigert polynomials were recently studied in [29] and [55], but 
in different asymptotic regimes: the degree goes to infinity while the weight 
is fixed. Therefore, the analysis of this section seems to yield new results for 
asymptotics of Stieltjes-Wigert polynomials. Nevertheless, the asymptotic 
analysis of the orthogonal polynomials and the ensemble (74) with varying 
weight (77) can be done in a very similar way to the analysis in [20, 21] using 
the Deift-Zhou steepest-descent method for related Riemann-Hilbert prob- 
lems (RHP's), which is now one of standard tools for asymptotic analysis of 
orthogonal polynomials. We note that [55] also used the Deift-Zhou method 
(for a different asymptotic regime) and our analysis has some overlap with 
the analysis of [55]. In this section, we present only a sketch of the analysis. 

It is a standard result in random matrix theory (see, e.g., [41, 52]) that 
(74) equals 



(80) 



det(l 



K fc W) 



where 



(81) K k (x,y) = yfw(x)w(y) 



Ik 



Pk+i{x)pk{y) - Pk(x)pk+i(y) 



x-y 
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is the Christoffel-Darboux kernel in which p n (x) = 7n£ n + • • • is the nth 
orthonormal polynomial with respect to w. Hence, 

E(F k z (Z))=det(l-K k H 1 ), 

(82) 

E(Gf(Z)) = det(l-K fc fla). 

Let Y(z) be the solution to the following Riemann-Hilbert problem: Y(x) 
is the 2 x 2-matrix-valued function on C \ i?+ satisfying 

• Y(z) is analytic for z G C \ M+, Y-j-(z) = lim e jo Y(z ± ie) is continuous for 
z G M + and Y(z) is bounded as z — ► 0; 

• for z G R+, 

(83) Y+(z) = Y_(z)(J ^ 

• Y(z)z~( fe+1 ) <T3 = (I + 0(z -1 )) uniformly as z -> 00 such that z G C \ 1+, 
where a 3 = (J _°) . 

There is a unique solution Y to this RHP and, in particular, the (11) 
and (21) entries of Y(z) are given by Yn(z) = 7^ 1 pfc + i(z) and Y2i(z) = 
—2ni'ykPk(z) [25]. Note that the existence of Y under the condition that 
Y(z) is bounded as z — ► (rather than, e.g., that Yi2(z) = 0(z~ 1 ), as in, 
say, [54]) is due to the fact that w(x) — > faster than any polynomial as 
x — > 0. Thus, the Christoffel-Darboux kernel can be written as, by using the 
fact that detY(z) = 1, 

(84) K k (x,y) = ^w(x)w(y) 2m{ 1 x _ y) (0 l)Y~\y)Y(x) f^. 

One of the main ingredient in analyzing the RHP for orthogonal poly- 
nomials asymptotically is the so-called equilibrium measure and the corre- 
sponding "^-function. " Let ip(x)dx be a measure on M + = supp(w) with 
total mass 



(85) / ip(x) dx = k + l. 
Define the "G-function" 

(86) G{z) = I log(z - x)i/)(x) dx, zGC\' 



where log represents the log function on the standard branch so that log u = 
log I it I + iarg(u), where | arg(n)| < n. It is customary to define ip to be the 
probability measure and define the g-function as in (86) [hence, G(z) = 
{k+ l)g(z)], but in this paper, we use the above convention since it simplifies 
some formulas below. Note that 

(87) G + {x) + G-(x) = 2 j\og\x-y\i>{y)dy, xGR + . 
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We look for G satisfying the following two conditions: there exists a constant 
£ such that 

• G+(x) + G-(x) + log(u;(a;)) — £ = for x G supp(V>), 

• G + (x) + gS{x) + \og{w{x)) -£<0 for x G K+ \ supp(^)- 

For such G, the measure ^ is called the equilibrium measure. 

Using the standard procedure to solve this variational problem (see, e.g., 
[19, 47]), one can compute the equilibrium measure for the weight (77). 

Lemma 4. For the weight (77), the support of the equilibrium measure 
is [a, b], where 



(2fc+l)/fc 2 _ ^/ e (4fc+2)/fc2 _ e 2/fc ; 
= e ( 2fc+1 )/ fc2 + ^/ e (4fc+2)/fe2 _ e 2/fc. 

77ie equilibrium measure is, for x G [a, b], 
(89) 

where R(z) = ((z — a)(z — b)) 1 ' 2 denotes the principal branch of the square 
root function and the simple closed contour C contains z and [a, b] , inside 
does not touch (— oo,0] and is oriented counterclockwise. A residue calcula- 
tion yields that 

, s , / » k 2 ( J(h — x)(x — a) \ . , . 

(90 ib(x) = arctan v v - A '- , x G fa, b . 

We remark that a and b are sometimes called the Mhaskar-Rakhmanov- 
Saff numbers. The above a and b are obtained in [55]: with a n and (3 n as in 
(2.2) and (2.3) of [55], we have 

(91) a=(e- 1 /( 2fc2 )a n )| fc ^ fc/2in=fc+1 , b = (e-V( 2 ^)/^)|^ fc/2>n=fc+1 . 

Given this ip, G(z) is defined as in (86) and I is defined as I = 2G(b) — 
log(w(b)) = 2G(a) — log(w(a)). The function h(z) in (89) is analytic in z G 
C \ (— oo,0]. A residue calculation yields that 

k 2 ^Jy/ab + z-R(zy 
r ab + z + R(z). 

where log denotes the principal branch of logarithm. For a computation 
below, we note that as k — > oo, 

(93) * 

v*=i + f k+ l + o {k -^. 



, \ , / x k 2 ( \fab + z - R(z)\ „ . . 
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We also remark that with x = 1 + ^7=, for w = 0(1), as k —> 00, at least 
formally, 

u 

(94) ip{x)dx V2-w 2 dw, we [-V2,V2], 

7T 

which is precisely Wigner's semicircle. This last calculation is not going to be 
used below, but it provides an intuitive reason as to why the ensemble (74) 
[and (72)] has the same asymptotics as the Gaussian unitary ensemble, not 
only locally, but also globally. 
Set 

(95) M(z) = e -(VW°3 Y ( z)e -G(z)v 3e (i/2)e* 3 

for z £ C \ 1R+ . Using the analyticity of G for z G M+ \ [a, b] and the varia- 
tional conditions, M(z) solves the following, equivalent, RHP: 

• ~M.(z) is analytic for z E C\R + , M.±(z) is continuous for z € M+ and M(z) 
is bounded as z — > 0; 

• for z G E+, M + (z) = M_(z)V M (z), where 

/ e G-(z)-G+{z) 1 \ 

(96) V M (z) = ( Q e G?+(*)-G_(*) ) > ^G(a,b), 

(97) V M (z) = L j J, ^R+\(a,b); 

• M(z) = I + 0(z~ 1 ) as z -> 00. 

The nonunit terms in the jump matrix can be expressed in a unifying 
way. Set 

(98) H(z) = G(z) + \\og{w(z))-\l, z€C\((-oo,0]U[a,b]). 
Noting the variational condition, we find that for z E (a, b), 

G+(z) -G-(z) = 2G + (z) +log{w{z)) -£ = 2H + (z) 

(99) 

= -(2GL(z) + log(u;(z)) -£) = -2£T_(z). 

Hence, G+(z) — G_(z) has an analytic continuation both above and the 
below the real axis. Therefore, the jump matrix Vm equals 

(100) V M (z) = [ Q e -2ff_wJ> ^G(a,b), 

(101) V M {z) = ( q ! J. zEM+\(a,b). 

Using the definition of G and Lemma 4, one can check that 

(102) H'{z) = \R(z)h(z). 
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We now scale the RHP for M so that the interval (a, b) becomes (—1, 1). 
In other words, instead of moving the interval as the support of the equilib- 
rium measure, we will fix the support. In that way, we can use the analysis 
of [20, 21] more directly. Define 



b — a b + a 



b + a 



(103) N(z)=M 
Set E = (-£tfi, oo) and set 

(104) H(z) = H 
The matrix N solves the following RHP: 

• N(z) is analytic for z G C \ E, N-t(z) is continuous for z G E and N(z) is 
bounded as z — > — ; 

• for ze S, N+(z) =N_(z)Vat(z), where 



(105) 



Vn(z) 



e ~-2H+(z) i 

e-^-W ^ ' 



z G 



■1,1] 



1 e 2#(*) 

1 



(106) Vat(z) 
• N(z) = I + 0(z _1 ) as z -> oo 



zGS\(-l,l); 



Note the factorization for z G (—1,1), 







1 

-2i3_ 0) 



1 

g -2H_( 2 ) x 



1 
-1 



1 

e -2H+(z) 1 



(107) 

where we use the fact that H+(z) + H_(z) = for z G (—1,1). Let Ej, j = 
0, 1, . . . , 4, and Qj, j = 1, . . . , 4, be the contours and open regions given in 
Figure 1. Contours are oriented from left to right. Define 

'N(z), zGfiiU0 4 , 
N ( z )(_ e -W) J), z£Q 2 , 

Then Q+(z) = Q_(z)Vq(z) for z in Eo, . . . , E4, where 

1' 



(108) 



Q(z) 



(109) 
(110) 
(111) 



Vq(z) 
Vq(z) 



-1 
1 

1 e 2^(^) 

1 



z G E , 
ZGE1UE2, 
zG E3UE4. 
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Fig. 1. Contours for N . 



The off-diagonal terms of Vq on Si U • • • U S4 converge to as the following 
lemma implies. 

Lemma 5. There exist 5q > and ko > such that for k>ko, 

(112) Re[H(x + iy)) > 2k\y\ \/l - x 2 for -1 < x < 1 and -S < y < S . 
For any 5 > 0, 

(113) H(x)<-k0 2 for -^-^ <x<-l- 5 and x> 1 + 5 

b — a 

when k>ko, and 

(114) lim / e 26{x) dx = 0. 
fc^oo J(s 3 us 4 )n{|^-i|><5}n{|2:+i|>5} 

Hence, Vq — ► Voo for a constant matrix Vqo defined as 

(115) V oo (z)=(^_ 1 J), zES , 

and Voo(z) = I for z G Si U • ■ • U S4, where the convergence Vq — > V m is in 
L°°(SoU---US4) and also in L 2 ((S U • • -U S 4 ) n {\z - 1| > 5}f}{\z + l\ > 5}) 
for an arbitrary, but fixed, 5 > 0. Let 

(116) «*> = U + l) ' 

where the branch cut is [—1, 1] and (5{z) ~ 1 as z — > +00 on the real line, 
and define 

(117) Q™(z)- l ( /3 + /3 " 1 -^-^)\ 

for z£C\ S . Then Q°°(z) is the solution to the RHP for the = Q^V^ 
and Q°°(,z) — > I as z —> 00. The convergence Vq — > V m is not uniform near 
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the points z = ±1, hence it is not true that Q(z) — ► Q°°(z) for all z and one 
therefore needs local parametrix for z in a neighborhood of ±1. 

Let ^(z) be the matrix- valued function constructed from the Airy func- 
tion and its derivatives, as defined in Proposition 7.3 of [20]. Let e > 0. For 
z G U r := {z : \z — 1| < e}, set 

(118) S r (z) = E(z)^((-|i?(z)) 2/3 )e-^ CT3 , 

where 

E(z ) = v ^ e (^(_ i i :^ 

(119) 




Note that E(z) is analytic in U r if e is chosen sufficiently small. The matrix 
Si(z) is defined in a similar way for z G U\ := {z : |z + 1| < e}. Define 

rQ°°(z), 2ec\f/ r u[/ius, 

(120) Q par (z) = I S r (z), zeU r \X, 

[s,(*)> z€Ui\S. 
From the basic theory of RHP, the estimate in Lemma 5 and the same 
argument as in [20], one can check that the jump matrix for Qp^Q converges 
to the identity in I? n L°°. Hence, 

(121) Q(z) = (I + 0(k~ 1 ))Q pSbr (z). 

This holds uniformly for z outside an open neighborhood of the contours 
£ U dU r U dUi. But a simple deformation argument implies that the result 
is extended to z on the contours (see [20]). Hence, by reversing the trans- 
formations Y -» M -» N -» Q [see (95), (103) and (108)], the asymptotics 
of Y(z) for all 2 G C are obtained. 

By substituting the asymptotics of Y into (84), edge and bulk scaling 
limits of the are obtained; see [15, 18, 21] for details. For xq such that 
Vk(xQ — 1) lies in a compact subset of (V~k(a. — 1), \/k(h — 1)), for all £, rj in 
a compact subset of R, 

(122) —!-—~Kk (x + -^,x + - S(£, rj) 
ip{x ) \ ip{x) ip{x )J 

in trace norm for £ , 77 G R, where 

x sin(-/r(£ — r?)) 

123 se,q)= 

7r(£ - 7?) 

Here, we may replace ip(xo) by Kfc(xo,xo). The error is 0(/c _1 ), uniformly 
for £,77 in a compact set. The convergence is also in trace norm in the 
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Hilbert space L 2 ((—rj,r])) for a fixed t] > 0. From (82), by taking xo 
the limit (30) in Proposition 2 is obtained. 
At the edge of the support of i/j(x), set 



(124) 

As k — > oo 
(125) 



5, 



-^b^a/i(6) 



2/3 



Bk 



K* 6 + 



6 + 



V2 



A(£,r?) 



in trace norm in the Hilbert space L 2 ((£,oo)) for a fixed £, where 

Ai(OAi'(r ? )-Ai'(e)Ai(, ? ) 



(126) 



is the Airy kernel. Hence, from (82), the limit (31) in Proposition 2 is ob- 
tained. 



4. Generalizations and discussions. We comment on three issues in this 
section: the case in which the moment generating function does not ex- 
ist, finite-dimensional distributions and the connection of this work to q- 
orthogonal polynomials. 

No moment generating function. In this paper, we have assumed the 
existence of the moment generating function for the random variable in- 
crements of nonintersecting random walks. This is simply to improve the 
estimates. For the case E|A A -'| 2+<5 < oo, 5 > 0, there is a version of the KMT 
theorem which provides analogous estimates to those used in Section 2. As 
one would expect for this case, must grow more rapidly in k. Another 
method for achieving results similar to those of this paper is to use Skorohod 
embedding in order to embed the nonintersecting random walks into Brow- 
nian motions. In order to achieve this, one must assume that E|A~/| 4 < oo. 

Finite- dimensional distributions. The results of this paper focus on the 
limiting distributions of nonintersecting random walks at the fixed time 
t = 1 . It is also interesting to consider finite-dimensional distributions of the 
process, that is, in the correct scaling t±, . .. ,t n E [1 — Ak" 1 / 3 , 1 + Ak^ 1 ^ 3 ], 
the finite-dimensional distributions of the fluctuations of the top random 
walk should converge to those of the Airy process. A similar, but differently 
scaled, result should also be true "in bulk"; see, for example, [1, 46, 53] 
and references therein concerning the Airy process and other processes from 
random matrix theory. The methods of Section 2 are certainly applicable 
to this problem, however, the convergence of the finite-dimensional distri- 
butions of the nonintersecting Brownian bridges to Airy /sine processes does 
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not follow immediately from the analysis of Section 3. However, one can use 
a different approach based on the method of Eynard and Mehta [14, 23]. In 
this approach, an inversion of a matrix is crucial. After the completion of 
the present paper, Widom communicated to the authors how to invert the 
matrix. Work in this direction will appear in a future paper. 

Stieltjes-Wigert weight and q-orthogonal polynomials. In Section 3, the 
Riemann-Hilbert problems for the orthogonal polynomials with respect to 
the Stieltjes-Wigert weight (79) was analyzed in the Plancherel-Rotach 
asymptotic regime. The analysis yields the asymptotics of the Stieltjes- 
Wigert polynomials in the entire complex plane. Since Stieltjes-Wigert poly- 
nomials are examples of g-polynomials, this result also yields an asymptotic 
result for certain q-polynomials. 

Acknowledgments. The authors would like to thank Percy Deift and 
Harold Widom for useful discussions. 
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